Temperature gradient driven lasing and stimulated cooling 
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A laser can be understood as thermodynamic engine converting heat to a coherent single mode 
field close to Carnot efficiency. From this perspective spectral shaping of the excitation light gen- 
erates a higher effective temperature on the pump than on the gain transition. Here, using a toy 
model of a quantum well structure with two suitably designed tunnel-coupled wells kept at different 
temperature, we study a laser operated on an actual spatial temperature gradient between pump 
and gain region. We predict gain and narrow band laser emission for a sufficient temperature gradi- 
ent and resonator quality. Lasing appears concurrent with amplified heat flow and points to a new 
form of stimulated solid state cooling. Such a mechanism could raise the operating temperature 
limit of quantum cascade lasers by substituting phonon emission driven injection, which generates 
intrinsic heat, by an extended model with phonon absorption steps. 

PACS numbers: 42.50.Pq 



Lasing is facilitated by implementing gain within an 
optical resonator. For stimulated emission based ampli- 
fication this requires an inverted population distribution 
between the amplifier quantum states pQ. In theoret- 
ical descriptions pumping to generate inversion can be 
well described by coupling reservoirs of different temper- 
ature to the pump, the injection and the lasing transi- 
tion [21 [3] . Practically this requires precise spectral de- 
sign of the pump radiation to overlap mainly with the 
pump transition with strong suppression on the lasing 
and injection transitions to yield inversion. Interestingly, 
such inversion can be consistently mimicked by simply 
assuming an effective negative temperature on the lasing 
transition [3]. At optical frequencies almost no thermal 
photons are present at room temperature and simple fil- 
tering of the pump light is sufficient. However, simple 
spectral shaping in the far infrared or even THz regime 
is practically hampered by ubiquitous thermal quanta, 
which create strong restrictions on the laser operating 
temperature [5]. 

In this paper we study possibilities to directly use spa- 
tial temperature gradients instead of spectral pump de- 
sign to generate inversion and facilitate gain and lasing. 
Although we mostly concentrate on the basic principles 
and mechanisms of such a setup, we also try to connect to 
a possible real world implementation based on designed 
quantum well structures. Here the optical properties and 
level structures in far infrared or even THz regime can be 
precisely engineered. While such a temperature gradient 
laser at first seems more of conceptual than of practical 
importance, the underlying principle can be easily gener- 
alized to a system with different pump steps. In partic- 
ular thermal or phonon absorption based injection steps 
to the gain transition lead to interesting modification of 
the thermal operation conditions of such a laser. Lasing 
can be tied to thermal absorption and heat flow to avoid 
unwanted heating and control heat flow [6]. 

The model we study is a simple structure comprised 



of two quantum wells coupled to a cavity. A sketch is 
shown in Fig.[l] The Hamiltonian describing a system of 
iV" two-well structures coupled to a single cavity mode is 
given by H = Ho + Hi with 
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The energy eigenstates are denoted by 
\uj) and \lj) with the eigenenergies tiujj 1 , 
and hidi and a' is the creation operator for a photon in the 
cavity mode with frequency u> m . The coupling between 
the u — ¥ I transition and the mode is described by the 
coupling rate g. Tunneling between the states \Ii j) and 
1/2,3) an d the states \I3.j) and is controlled by the 

tunneling matrix elements J\ and J2. For the forthcom- 
ing analysis we define all frequencies with respect to uji 
and switch to a frame rotating with w m introducing the 
detuning A M = cj„ — uj m . The dynamics of the reduced 
density matrix is described by the master equation 
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where the Liouvillian C describes dissipative processes, 
including the coupling to the environment at finite tem- 
perature, see supplementary material. 

The decay rate of the cavity is denoted by k, while the 
spontaneous emission rates of the remaining transitions 
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FIG. 1. Sketch of the two-well structure and the energy levels. 
The wave functions are assumed to be well localized in each 
well while the coupling is described by the tunneling constants 
Ji and J 2- 



are denoted by 71, . . . , 74. The thermal photon number 
at the transition frequencies oj\ = w u —wi, u}% = wj 3 —uj 2 , 
cj 3 = ujj i — tjj u and W4 = Wj — Wjj at temperature T is 
denoted by n t (T) (i = 1, . . . , 4). 

We assume that we have a temperature gradient across 
a sample containing N two-well structures and that T\ > 
T2 in each structure. As T\ ^ we will find population in 
the state I/3). Via the tunnel-coupling J\ the state 
gets populated as well. Due to the lower temperature T 2 
of the second well the population is preferably transferred 
to the state |it). With a certain probability a photon is 
created in the cavity and leaves it. The cycle ends with 
the tunneling from state to li^)- 

Using Eq. ([3]) we calculate operator expectation value 
equations for the populations in each state and the quan- 
tities they couple to [7J. The equations are truncated to 
include only second order correlation functions. The pop- 
ulation of the state \v,j) is described by 
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where we omitted the temperature dependence of the 
rii and used the assumption of total phase invariancc 
((a) = (\uj) = 0) to simplify the expansion of third 
order expectation values [S] . This leads to the next equa- 
tion 
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(a\uj) = iA u (a\uj) 

-i fl ([l + (at a >] (|^)( % |)-(ata><|^)fe|)) 
- (/e + 71 (2ni + 1) +73^3 +74 (ra 4 + 1)) (a\uj) 
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which describes the rate at which excitations are trans- 
ferred between the cavity and the structure. The num- 
ber of photons follows ^ (et^a) = — 2«^a^a) + 2/cni — 
igN (Ja^ \lj) (uj\^ — (a\uj) In the last line of 

Eq. aty we see that cooperative effects between the struc- 
tures emerge as a result of the common interaction with 
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FIG. 2. (Color online) (a) Numerical solution for the steady 
state number of photons (a'a) st in the cavity as a function of 
Ti and T2. The number of structures was chosen to N = 10 7 . 
The dashed line indicates Ti jt h necessary to obtain cooper- 
ative behavior (valid for T2 small enough to justify setting 
ni = nz = n 4 = 0). (b) (tr a) as a function of Ti and N 
(T 2 = 0.1 K). The analytical solution gives N th (Ti) (dashed 
line) and reliably predicts the onset of cooperative behavior 
and hence a significant increase in (a/a) . 



the cavity. We find 
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In addition we calculate equations for 

(i = 1,2,3,4), (I 2J \) and (|/ 3j ) (hjl)- 

Approximate analytical solutions can be deduced un- 
der certain assumptions and give insight into the main 
mechanisms [5]. We assume that T 2 is small enough to 
justify setting ri\ = 713 = 714 = and that A„ = 0. By 
keeping only the decay term proportional to k and the 
last line in Eq. ([5]), where we approximate N — 1 sw N, 
we obtain the steady-state 

(a\u j ){l j \) st = - 1 -^{\l j ){u j \\u k ){l k \) st . (7) 

Under the assumption of negligible T2 and by using 
Eq. ([7]) it is possible to calculate approximate steady 
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state solutions for all of the involved quantities. The 
steady state of the correlation between the structures 
( u j \ \ u k) ('fe|) s t ^ s °f particular interest. It shows for 
which parameters the formation of cooperative behav- 
ior with the mode is possible [5]. For a given number 
of structures N we can calculate the critical tempera- 
ture Ti )t h that is necessary to meet the condition for 
cooperative behavior (uj\ \uk) (h\) st > 0). At the 
same time, for a given temperature T\ and the resulting 
number n 2 (Ti) a certain threshold number JV t h of struc- 
tures is necessary for cooperative behavior. Details can 
be found in the supplementary material. 

The full set of equations is also integrated numeri- 
cally to determine the exact steady states and to com- 
pare them to the approximate solutions. The parameters 
K = 2tt x 5 x 10 5 Hz (i.e. FWHM linewidth of 1 MHz), 
g = («/3) x 1CT 3 , A u = ,71 = 2tt x 4 x 10 2 Hz, 72 = 
73 = 74 = 7i x 10 2 , ui = 2?r x 7.9 THz, w 2 = 2vr x 9 THz, 
w 3 = 2n x 0.1 THz, w 4 — 2n x 1 THz remain fixed for the 
forthcoming calculations, while N, T\ and T 2 are varied. 
The number of photons as a function of T\ and T 2 for 
N = 10 7 is shown in Fig. [2? a). The occupation of the 
cavity mode strongly depends on the temperature gra- 
dient between the wells. The approximative analytical 
treatment for small T 2 gives the threshold temperature 
T\ th that signals the onset of cooperative behavior. For 
T 2 > 1 K the approximative solutions break down. In 
Fig. p[b) we show (aJa) as a function of T\ and N 
for fixed T 2 = 0.1 K. Here, the approximative analyti- 
cal solutions give a reliable result for which parameters 
to expect a high number of photons in the cavity. The 
numerical solutions show that the number of photons in 
the cavity increases rapidly as soon as we exceed N t h for 
a given value of Ti . 

Assume that we have a layer containing N structures 
which initially is hot. If one side of the layer is cooled 
by some means, the resulting temperature gradient will 
lead to a population of the cavity mode. The rate at 
which photons are lost from the cavity is 2k (aJa). This 
suggests that the cavity can serve as a channel for cool- 
ing. Without the second quantum well, or equivalently 
with J\ = J 2 = 0, we find the ratio between the upper 
level -Z3 and the lower energy level J 2 in the first well to 
be 712/(1 + n 2j- In the presence of the second well, and 
in the regime where we find a considerable number of 
photons in the cavity, this ratio is reduced. Figure [3^ a) 
shows how the ratio (h,j\) I (|^2j) (^2,3 1) changes 

with increasing T\ while T 2 = 1 K remains fixed. With 
increasing 7\ the relative occupation of the upper level 
|/3j) increases but remains below the result for the un- 
coupled case. The difference starts to show at Ti ps 75 K 
which is the critical temperature above which we obtain a 
significant increase in (a^a) st , as can be seen in Fig. 21 a). 
The same comparison can be made for fixed T\ = 400 K 
(leaving n 2 /(l + n 2 ) fixed) and varying T 2 . While T 2 
remains low enough to ensure a significant number of 
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FIG. 3. Occupation ratio in the first quantum well for the 
coupled and uncoupled case versus Ti and T2. In (a) we 
choose T2 = 1 K and vary Ti. The ratio between the pop- 
ulation of I3J and I2,] (solid line) increases with Ti but re- 
mains below the dashed line which gives the ratio in the un- 
coupled case. In (b) the temperature of the hot part Ti is 
fixed to 400 K. For T2 < 9 K we again see the reduction of 

(l^)(/ 8J |)/(|/2,i></ 2 ,i|>. 



photons in the cavity, and hence a considerable dissipa- 
tion via the cavity mirrors, the occupation ratio in the 
first well is reduced. This corresponds to an effectively 
lower temperature of the first well. 

We can obtain a large number of photons in the cavity 
for N > A/th (Ti) and an appropriate ratio between Ti 
and T 2 . The temperature gradient can therefore be seen 
as an effective pump. To characterize the emission we 
calculate the spectrum of the cavity output field and the 
linewidth. Using the quantum regression theorem we ob- 
tain the differential equation for the first order correlation 
function (a* (t + r) a (i)) and (\uj) (t + r) a (t)), sec 
supplementary material. To find the spectrum we calcu- 
late the Laplace transform of the coupled equations [TD] . 
The steady state values of the involved quantities are ob- 
tained numerically. In Fig. [4] the linewidth (FWHM) of 
the cavity output spectrum is depicted. With the onset 
of cooperative behavior the linewidth is reduced dramat- 
ically. Similar behavior is found for fixed T 2 and varying 
Ti and N. 

The two-well system discussed so far clearly exhibits 
the essential physical mechanism of heat to coherent light 
conversion in its purest form. In the following we show 
that the principle can be extended to multi step exci- 
tation processes, where several heat absorbing steps are 
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FIG. 4. (Color online) Linewidth of the cavity output spec- 
trum for varying Ti and T2 with N = 10 7 (plotted on a log- 
arithmic scale). The linewidth (FWHM) of the empty cavity 
2^ = 1 MHz corresponds to the dark red color (upper bound- 
ary of the scale). The linewidth A/ falls slightly below that 
limit with increasing Ti until it finally drops off dramatically 
as the cooperative effects emerge. For small T2 the onset of 
cooperative behavior is marked by Ti lt h- 



combined to lead to operation at higher emission fre- 
quencies. Mathematically, this can be easily facilitated 
by incorporating an absorptive step in the second well 
by choosing u/ u > uj±, without needing a more complex 
model. In the dissipative dynamics we simply have to 
exchange the rates 73713 (T2) «-> 73 (n 3 (T 2 ) + 1). Phys- 
ically, this means that in each cycle in this system the 
absorption of an extra thermal photon or phonon with 
frequency U)' u — uji i is necessary. This implies that we find 
a minimal temperature below which not enough phonons 
are present and lasing is impossible. In general one could 
of course even think of a combination of thermally as- 
sisted and externally driven excitation steps. We will re- 
frain from this here to keep the complexity of the model 
limited. 

Let us look at the results now. The temperature 
range for which (a 1 a) becomes significant is depicted 
in Fig. [5]ja). The dashed line marks the critical tem- 
perature necessary for cooperative behavior, which is ob- 
tained analytically with similar approximations as de- 
scribed above. The linewidth of the emission is shown in 
Fig. [5jb) . This concept could be incorporated as an extra 
step in conventional quantum well lasers, which generate 
intrinsic heat during operation, to extend their range of 
operation as the absorption of phonons would counteract 
the heating process. 

Our theoretical calculations exhibit the possibility to 
use a spatial temperature gradient in a suitable active 
medium to induce optical gain and lasing. While not very 
practicable in the visible regime, a semiconductor het- 
erostructure can be envisaged to generate narrow band- 




600 800 

(b) 



1000 1200 1400 
Ti[K] 




< 



aC 
O 

2 " 



600 800 1000 1200 1400 
Ti[K] 



FIG. 5. (Color online) (a) Number of photons in the mode 
versus Ti and Tb- For T2 > T^th we find a sudden increase 
in (a'a) , and in (b) a corresponding drop in the linewidth 
of the cavity emission. The parameters remain unchanged 
except for 73 =71 = 2n x 4 x 10 2 and u' u = ui u + 2tt x 0.2 THz. 



width emission through a resonant cavity mode in the 
THz regime, where the coherent emission of photons is 
sustained by the spatial temperature gradient. The ef- 
fective extraction of thermal energy from the system in 
parallel provides for stimulated cooling of the first, hot 
quantum well. To use a similar idea at higher frequen- 
cies without too high temperatures, a multi step excita- 
tion process can be envisaged, where some steps could 
also invoke optical pumping or current injection. In gen- 
eral, introducing an absorptive step well before carrier 
injection to the upper laser level relies on absorption of 
thermal quanta and therefore counteracts heating during 
operation. This leads to a higher operating temperature 
limit. While we here only worked with an oversimplified 
generic model, the physical principles behind should stay 
valid for more realistic and thus more complex descrip- 
tions. In addition, transverse spatial temperature gradi- 
ents or even spatially separated structures connected by 
conducting wires might provide for alternative configu- 
rations at greater length scales. 
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SUPPLEMENTARY MATERIAL 
Liouvillian 

The Liouvillian that describes dissipative processes, in- 
cluding the coupling of the cavity and the structures to 
the environment at finite temperature, reads 



(8) 



Approximate Analytical Solutions 

For T 2 small enough we can neglect rii, n 3 and 
714. Inserting the steadystate solutions for (a\uj) (lj\) st , 
(\uj) (uj\) st and (\lj) (lj\) st into Eq. ^ we obtain 



C [p] = K (rii (T 2 ) + 1) (2apa) — a)ap — pa) a) + tm\ (T 2 ) (2a) pa — aa) p — paa)) 

+ 71 («i (T 2 ) + 1) 53 (2 (uj\ p \uj) - \u 3 ) ( Uj \p - p \uj) (uj\) 
3 

+ 7i (nx (T a )) ]T (2 \ Uj ) (lj\p ( Uj \ - \l 3 ) (l 3 \p - p \l,) {l 3 \) 
3 

+ 72 (na (TO + 1)£ (2 |7 2ii ) (I 3d \ p \I 3J ) (I 2 . 3 \ - \I 3<j ) (I 3J \p - p \I 3<j ) (I 3J |) 
3 

+ 72 (na (TO) ]T (2 |I 3J ) (I 2ji | p |J 2)i ) (I 3J | - \h,j) (hj\p - P \h,j) (h,j\) 
3 

+ 73 (n 3 (T 2 ) + 1) (2 \ Uj ) (h d \p \I id ) (u 3 \ - \h d ) (h i3 \ p-p \I A>S ) (I 4 , 3 \) 
3 

+ 73 (n 3 (T 2 ) (T 2 )) ^ (2 |J 4J ) (u 3 \ p \u 3 ) (I 4ji | - \u 3 ) (u 3 \ p - p \u 3 ) (u 3 \) 
3 

+ 74 (n 4 (T 2 ) + 1) J2 ( 2 \hj) (hi P \h) (hj\ - Up - P \h) M) 
3 

+ 74 (n 4 (T 2 )) ( 2 \h) (h,j\p \h,j) ~ \hj) (hjp - P \hj) (hj\) 
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g 4 N 2 | J 2 | 2 ( Uj \ \u k ) (l k \) st (2 + 3n 2 ) 7273 ( 72 + n 2l2 + 73 ) 74 + | Ji| 2 (g A N 2 ( Uj \ \u k ) (l k \) st n^hs!* 

+ |^| 2 (2g 4 N 2 (uj\ \u k ) (l k \) st ((72 + 73)74 + "-272(73 + 274)) +g 2 ^27273(71 - 74 + "-27273(74 - 7i) k2 )) 

/( 7ik 2 ( I J 2 | 2 (2 + 3n 2 ) 7273 (7 2 + n 272 + 73 ) 74 + | Ji| 2 (" 272 7374 + 2 | J 2 | 2 ((72 + 73)74 + "272(73 + 274)))) 



(9) 



for the quantity (\lj) (uj \ \u k ) (^fc|) st describing coopera- closed by the rectangular brackets. To find the number 
tive behavior between the structures. The stable solu- of structures necessary for reaching the onset of cooper- 
tion is obtained from the part of Eq. ^ which is en- ative behavior iV t h (Ti) for a given temperature Ti we 

solve (uj\ \u k ) (l k \) st = to find 



N t h (Ti) = (71 + 74 )k \J 2 \ (2 + 3n 2 ) 7 i 7273 (7 2 + n 272 + 73 ) 74 



+ \ J i\ 2 ("271727374 + I (27i(72 + 73)74 + "272(7173 + 4 7 i 7 4 +7:374))) 



/[ g 2 "27273(74 - 7i) ki| 2 l^| 2 I • (10) 



The result is plotted in Fig. [2jb) (dashed line) . Similar 
calculations give the critical temperature Ti jt h for given 
N, see Fig. Ha). 



Calculation of the spectrum 

Using the quantum regression theorem we obtain 
(a* (t + r)a (£)) which couples to (\uj) \ (t + t) a {£)). 
For t — > 00 using stationarity we can write 



d / (at(r)a(O)) 
dr \(\u 3 )(l 3 \(r)a(0)) 



— K 



-igN 



-i5((K(%l)) st -(l^(^l» st ) iA„-r; V(K)fe|(r)a(0)) 



(at (r)o(0)) 



(11) 



with r = 7 3 " 3 (T 2 ) + ji (2ni (T 2 ) + 1) + -74 (71,4 (T 2 ) + 1). which allows to calculate the spectrum |10) . 
The Laplace transformation yields 



<ato) (a) = 



(g f g) st (s - iA u + r) - igN (\ Uj ) (l 3 \ a) st 
(s + «)(*- iA u +r) - g 2 N ({\ Uj ) ( Uj \) st - (\l 3 ) &I) J 

(12) 



